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Abstract 



a.- 

Conventional and hybrid light quark baryons are constructed in the non-relativistic flux-tube 



model of Isgur and Paton, which is motivated by lattice QCD. The motion of the flux tube with 
the three quark positions fixed, except for center of mass corrections, is discussed. It is shown that 



> 

| the problem can be reduced to the independent motion of the junction and the strings connecting 

o : 

t^. ■ the junction to the quarks. The important role played by quark-exchange symmetry in constraining 

o : p 

CNJ . the flavor structure of (hybrid) baryons is emphasized. The flavor, quark spin 5 and J of the seven 

O 1 

low-lying hybrid baryons are found to be (N, A) 2S+1 J P = iV 2 ± + , iV 2 § + , A 4 i + , A 4 | + , A 4 f + , 
where the iV 2 i + and iV 2 ! + states are doublets. The motion of the three quarks in an adiabatic 

o : 

potential derived from the flux-tube dynamics is considered. A mass of 1870 ± 100 MeV for the 



lightest nucleon hybrids is found by employing a variational method. 
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I. INTRODUCTION 



Since all possible good (J p ) quantum numbers of baryons can be described by conven- 
tional excitations of three quarks, the description of hybrid baryons, defined as bound states 
containing explicit excitations of the gluon fields of QCD, is necessarily model dependent. 
Nevertheless, any model of QCD bound states which allows the gluon fields to be dynamical 
degrees of freedom, as opposed to simply generating a potential (or surface) in which the 
quarks move, will have additional states involving excitations of those degrees of freedom. 
A description of the spectrum of hybrid baryons, the degree of mixing between them and 
conventional qqq excitations, and their strong decays will therefore be necessary in order to 
describe the results of scattering experiments which involve excited baryons. For example, 
such experiments make up the excited baryon resonance (N*) program at Jefferson Lab- 
oratory, where many excited states of baryons are produced electromagnetically. Hybrid 
baryons must play a role in such experiments. In principle their presence can be detected by 
finding more states than predicted in a particular partial wave by conventional qqq models. 
Doing so will require careful multi-channel analysis of reactions involving many different 
initial and final states 111. Another possibility is that such states will have characteristic 
electromagnetic production amplitudes |2J. If hybrid baryons obey similar decay selection 
rules to hybrid mesons |3[, they may be distinguishable based on their strong decays. This 
work concentrates on a determination of their masses and quantum numbers, in an approach 
where the physics of the confining interaction defines the relevant gluonic degrees of freedom. 

Hybrid baryons have been examined using QCD sum rules |^| and in the large number of 
colors (large N c ) limit of QCD One approach to modeling the structure of hybrid baryons 
(not taken here) is to view them as bound states of three quarks and a 'constituent' gluon. 
Hybrid baryons have been constructed in the MIT bag model JfJ by combining a constituent 
gluon in the lowest energy transverse electric mode with three quarks in a color-octet state, 
to form a color-singlet state. With the assumption that the quarks are in an S -wave spatial 
ground state, and considering the mixed exchange symmetry of the octet color wavefunctions 
of the quarks, bag-model constructions show that adding a J p = 1 + gluon to three light 
quarks with total quark spin 1/2 yields both N (I = |) and A (/ = |) hybrids with 
J p = | + and | + . Quark spin 3/2 hybrids are N states with J p = | + , | + , and | + . Energies 
are estimated using the usual bag Hamiltonian plus gluon kinetic energy, additional color- 
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Coulomb energy, and one-gluon exchange plus gluon-Compton 0(a s ) corrections. Mixings 
between q 3 and q 3 g states from gluon radiation are evaluated. If the gluon self energy is 
included, the lightest N hybrid state has J p = \ + and a mass between that of the Roper 



resonance and the next observed J p = \ + state, iV(1710). A second J p = \ + N hybrid and 
a J p = | + N hybrid are expected to be 250 MeV heavier, with the A hybrid states heavier 
still. A similar mass estimate of about 1500 MeV for the lightest hybrid is attained in the 
QCD sum rules calculation of Ref. |4|. 

For this reason there has been considerable interest in the presence or absence of light 
hybrid states in the Pn and other positive-parity partial waves in nN scattering. Interest- 
ingly, quark potential models which assume a q 3 structure for the Roper resonance j?l 0| 
predict an energy which is roughly 100 MeV too high, and the same is true of the A(1600), 
the lightest radial recurrence of the ground state J p = | + A(1232). Furthermore, models 
of the electromagnetic couplings of baryons have difficulty accommodating the substantial 
Rope, r _e photoconpling extracted from pion production data Q. Evidence 
for two resonances near 1440 MeV in the Pn partial wave in irN scattering was cited [10J, 
which would indicate the presence of more states in this energy region than required by 
the q 3 model, but this has been interpreted as due to complications in the structure of the 
Pn partial 
lations 



wave in this region, and not an additional qqq excitation |ljj. Recent calcu- 
li of Ntt — > Ntttt reaction observables incorporating baryon-meson dynamics are 
able to describe this reaction in the Roper resonance region in the absence of a qqq excita- 
tion, and find a dynamically-generated pole at the mass of the Roper resonance. Given this 
complicated structure, it is perhaps not surprising that there are difficulties in describing 
the photocouplings of this state within a simple three-quark picture. 

The motivation of this work is to build a model consistent with predictions from QCD 
lattice gauge theory, based on the Isgur-Paton non- relativist ic flux-tube model This 
model is motivated from the strong coupling limit of the Hamiltonian lattice gauge theory 
formulation of QCD (HLGT). This strong coupling limit predicts linear confinement in 
mesons proportional to the expectation of the Casimir operator for color charges, which has 



been verified in lattice QCD 14j. In conventional baryons, in the limit of heavy quarks, 
the static confining potential has been shown in lattice calculations to be consistent 
with that given by a minimum-length configuration of flux tubes meeting in a Y-shaped 
configuration at a junction (see Fig. QJ, and not consistent with two-body confinement, 



3 



where a triangle of tubes would connect the quarks in a A configuration [if]]. It is possible 
to experimentally examine this configuration in studies of baryon production in the central 
rapidity region of ultrarelativistic nucleon and nuclear collisions |17| . 

This structure of the glue, where the gluon degrees of freedom condense into flux tubes, 
is very different from the constituent-gluon picture of the bag model and large- iV c construc- 
tions. Substantial progress has been made in recent years in understanding conventional 
baryons by studying the large N c limit of QCD. However, the large N c limit does not nec- 
essarily provide model-independent results on hybrid baryons. Hybrid baryons in the large 
N c limit consist of a single gluon and N c quarks [5j. Even in the case of physical interest, 
N c = 3, this does not correspond to the description of hybrid baryons presented here, where 
is is argued that the dynamics relevant to the structure of hybrids are that of confinement, 
where numerous gluons have collectively condensed into flux tubes. Since the color structure 
of a hybrid baryon determines (through the Pauli principle) its flavor structure, the N c = 3 
case gives rise to states with different quantum numbers than in our approach. It has been 
shown by Swanson and Szczepaniak |l8| that a constituent-gluon model is not able to re- 
produce lattice QCD data for hybrid-meson potentials at large interquark seperations. 
In addition, the flux-tube model hybrid-meson potential is consistent at large interquark 
separations with that evaluated from lattice QCD 13, 2(3 ] . 



Hybrid baryons are constructed here in the adiabatic approximation, where the quarks 
do not move in response to the motion of the glue, apart from moving with fixed interquark 
distances in order to maintain the center-of-mass position. The effect of the motion of the 
glue in hybrid baryons (and the zero-point motion of the glue in conventional baryons) is to 
generate a confining potential in which the quarks are allowed to move. This differs from 
that found from multiplying the sum of the lengths of the strings ("triads") connecting the 
quarks to the junction by the string tension. The adiabatic approximation is exact only in 
the heavy quark limit, although the success of quark model phenomenology of conventional 
mesons and baryons implies that there is a close relation between heavy-quark and light- 
quark physics. A modified adiabatic approximation is employed, which can be shown to give 
exact energies and wave functions for specific dynamics even for light quarks |2l( . Moreover, 
a modified adiabatic approximation has been shown to be good for properties of light quark 
mesons in the flux-tube model [23j . 

The model is motivated from the strong coupling limit of HLGT, where there are "flux 
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FIG. 1: A possible configuration of flux lines in a baryon on the lattice. 



lines" which play the role of glue. In the adiabatic approximation operators which make 
the quarks move are neglected. The plaquette operator corrects the strong coupling limit, 
and induces motion of the "flux lines" between the quarks and the junction perpendicular 
to their rest positions. The flux lines are modeled by equally spaced "beads" of identical 
mass, so that the energy of each flux line is proportional to the number of beads, and hence 
its length. The spacing of the beads along the rest positions of the flux lines can be thought 
of as a finite lattice spacing, and the beads are allowed to move perpendicular to their rest 
position. The beads are attracted to each other by a linear potential, and the resulting 
discretized flux lines vibrate in various modes. Global color invariance requires that the 
three flux lines emanating from the quarks meet at a junction, which is also modeled by 
a bead. It is shown in HLGT that a single plaquette operator can move the junction and 
retain the Y-shaped string with the links in their ground state (see Fig. Eland Appendix |A"|). 
so that the junction may have a similar mass to the beads. However, for generality we allow 
the junction to have a different mass associated with it than that of the other beads. 

The final picture of both conventional and hybrid baryons is that of three quarks, con- 
nected via a line of beads to the junction in a Y-shaped configuration. The potential between 
neighboring beads is linear. The adiabatic approximation is used, so that the string is as- 
sumed to adjust its state quickly in response to motion of the quarks, thus generating a 
potential in which the quarks move, in both conventional and hybrid baryons. The motion 
of the quarks in these potentials is then solved for variationally. A brief outline of this 
approach is given in Ref. 2^]. The purpose of the present paper is to describe the model 
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FIG. 2: Flux lines in a baryon on the lattice, illustrating the application of a pair of plaquette 
operators, the upper operator moving one flux line perpendicular to its starting position, and the 
lower attempting to move the junction. 



and the calculation of hybrid baryon masses in much more detail, and to put this work in 
the context of recent advances in lattice gauge theory. 

The rest of the paper is organized as follows. Section |H] describes the dynamics of the flux 
tubes in various quark configurations, with analytic solutions in special cases. Section IHII 
discusses the quark-label exchange symmetry, parity, and chirality of the flux configuration. 
In Section ITVl the orbital angular momentum and color of the flux, and the combined quark 
and flux wave function are constructed. The following Section IS describes the potential 
in which the quarks move, which includes the energy of the flux. Numerical estimates of 
the masses of hybrid baryons are given in Section IV1I In Sections IV 111 and IV 1111 further 
discussions and conclusions are given. 



II. FLUX DYNAMICS 



Denote by #132 the angle between the line from quark 1 to 3 and that from quark 2 to 3. 
If #123, #132, and #213 are all smaller than 120°, the flux is in its equilibrium (lowest-energy) 
configuration when the junction is located such that there are angles of 120° between each 
of the triads that connect each of the quarks to the junction, and the beads all lie on the 
triads. In this lowest energy configuration, the string lies in the plane defined by the three 
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TABLE I: Notations frequently used in the main text. 



b string tension 

8 j angle between triads suspended at the equilibrium junction position 

i = 1,2,3 quark or triad label 

Ni number of beads on triad i 

n = l,Ni counts the beads on triad i from quark i to the junction 

m = 1, Ni counts the modes of the triad i 

Mi mass of quark i 

m& mass of the beads 

raj mass of the junction 

li distance from the equilibrium junction position to quark i 

li direction from the equilibrium junction position to quark i 

ri position of quark i 

ql amplitude of mode m on triad i in the QQQ plane, but perpendicular to h 

Q±m amplitude of mode m on triad i perpendicular to the QQQ plane 

r = (x, y, z) cartesian coordinates of the junction 

i) z junction oscillation directions parallel and perpendicular to the QQQ plane, respectively 

u±, lu z junction oscillation frequencies parallel and perpendicular to the QQQ plane, respectively 




FIG. 3: Flux configuration when none of the angles in the triangle joining the quarks are larger 
than 120°. 

quarks, denoted by the QQQ plane. The angle between the line from any two quarks to 
the junction equilibrium position r e q is 120°, which is denoted by 9j = 120° (see Fig. Eland 
Table H}. In terms of the lengths U of the lines from the i-th quark to the junction, and the 
quark positions r^, the equilibrium junction position is 

^ hh + V3 + hh 

The position vectors of the quarks relative to the junction equilibrium position are li = 
?i ~ r e q- 
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FIG. 4: Flux configuration when one of the angles in the triangle joining the quarks, here #132, is 
larger than 120°. 
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FIG. 5: Coordinate system used to describe motion of the flux in the QQQ plane, for configurations 
with (a) all of #123, #132 , and #213 less than 120° and (b) with 0j = #132 larger than 120°. 

If one of 0i23, #132, and #213 is larger than or equal to 120°, the equilibrium configuration 
of the flux is not this Y-shaped configuration. If 9^ > 120°, the lowest energy configuration 
has the junction at the position of quark j (see Fig. UJ). This situation is denoted by 
9 j > 120°. In what follows the case where 8 132 > 120° is analyzed, but the formulae for the 
other cases follow by the appropriate label exchange. 

An axis system is chosen as indicated in Figure El This defines normalized x and y, which 
can also be written in the 9j = 120° case as 

l± — I2 a ^1 + ^2 — 2Z3 



V3 



(2) 



where the are unit vectors along the triads, and h + 1 2 + h = 0, so that y equals l 3 and x 
is perpendicular to Z3. The third triad lies on the positive y-axis and the other two triads are 
120° on either side of the y-axis: triad one on the left-hand and triad two on the right-hand 
side. It is assumed that the sum of the masses of the beads and the junction is the energy 
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FIG. 6: Motion of the flux configuration decomposed into "junction motion" and "bead motion". 



of the flux configuration in its equilibrium position, so that 



3 3 



which implies 



i=l 



i=l 



V d N- + ^ 



(3) 



(4) 



where iVj is the number of beads on triad i. Note from the above that the bead mass m& is 
determined by the string tension b, the triad lengths k, and the ratio mj/rrib, and so should 
not be regarded as an independent parameter. 



A. Hamiltonian for 6j = 120° case 

The flux configuration is made dynamical by allowing the junction and the beads to 
vibrate with respect to their equilibrium configuration. There are two important motions 
which are expected to have physical significance: (1) the motion of the junction perpendicu- 
lar to and within the plane relative to its rest position, denoted "junction motion", and (2) 
the motion of the beads in the two directions perpendicular to the line connecting the quark 
to the junction, called "bead motion" , as illustrated in Fig. 13 The bead motion coordinates 
are not their positions, but the oscillating- wave amplitudes (defined in Appendix [Bj of the 
beads relative to their rest positions on the triads. It is important for what follows that the 
bead position coordinates are defined relative to their rest positions on the triads between 
the quarks and junction, which have followed the junction motion (see Fig. Ej). The Hamil- 
tonian is written in terms of the junction and bead motion coordinates. In what follows the 
small-oscillations approximation is used, where the beads and junction remain close to their 
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FIG. 7: Coordinates used to describe the motion of the beads in the presence of a displaced 
junction. 

positions in the equilibrium configuration. This approximation is used to motivate the basis 
of the subsequent numerical treatment, which is that it is a reasonable approximation to 
treat the flux motion as that of the junction, with an effective mass which depends on the 
equilibrium lengths of the triads, among other quantities. In the numerical treatment that 
follows, the restriction to small oscillations is removed. 

As the string moves, the quarks are allowed to move with fixed positions relative to each 
other, in order to keep the center of mass fixed. This is called the "redefined adiabatic" 
approximation. By working in this approximation, some but not all non-adiabatic effects 
are incorporated. 

The flux Hamiltonian for the 9j = 120° case in the redefined adiabatic approximation is 
(see Appendix iBj) 



tfflux = o^eff^ + V! 



sm.osc. 



i=l m=l 

3 N, 



h 3 Ni 1 

2 5>E2 ( ^ )2 Mrn? + (Vim?] 
i=l m=l 



mi 



=1 m=l 
1 



0, 



m b 

o ^-^3 



ELi Nk + 2 



2 ELiibh + Mj 



(4 m) 



3 Ni 



2 n=i(«*+K 



^ ^ a n a m{tf\\n<l\\rri + QlnQLr 



3 N, Nj 

ij^j n=l m=l 
3 N, 



O 



i=l n^m 



=1 m=l 



(at) 2 - 



K J 2 + telJ 2 ] 



(5) 
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The first term in Eq. El is the kinetic energy of the junction, with an effective mass of 



M. 



eff = m b 



E 



N^ + l) + mj_ 



m b - 



i 6(JV f + l) 



m 6 



4(5Xi(Wi + Mi)) 



(6) 



where the last term in Eq. El arises from the center-of-mass correction, and the first from the 
trivial motion of the beads which accompanies motion of the junction. The second term in 
Eq. El is the potential energy of the junction in the small-oscillations approximation, given 
in terms of the coordinates r = (x, y, z) defined in Fig. El by 
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yjm.OBC. =h^k 



i=l 



b 

+ 2 



x~ 



1 1 1 



-v 



l l 

h + T 2 



v/3 (\ 1\ 2 fl 1 1 

-T xy lk-T 2 +z [k + T 2 + T 3 



•(7) 



The third and fourth terms in Eq. El are the kinetic and potential energies of the beads, 
respectively, written in terms of the effective mass of the beads, including a center-of-mass 
correction, 



and the frequency of the m-th normal mode on the i-th triad, 

2(N t + 1) 



■ sin ■ 



(8) 



(9) 



k ""2(^ + 1) 

The fifth and sixth terms in Eq. El represent interactions between the junction and the 
beads, and interactions among the beads, respectively. In these terms the vectors e l are 
defined to be perpendicular to Zj, so that 



2' 2 ' 



1 -2^o 

2' 2 ' 



1,0,0), e *-et = -±. (10) 



In terms of the amplitudes q % m of the m-th normal mode on the i-th triad, where displace- 
ments of the beads along the triads are not allowed, the coordinates q^ m = e l ■ q l m are 
the projections in the QQQ plane, and the coordinates q l ±m are the projections of these 
amplitudes out of the QQQ plane. Finally, the (3 l m and a l m are the sums 

vA V ■ mpn ( — \\ m + l mn 



p=l 



2(iV, + 1) 
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and 



=E 

p=i 



sin 



mpir 







cot 



even 
odd. 



(12) 



2(iVi+l) 

The above demonstrates that the Hamiltonian can be seperated into three parts. The 
first (the first two terms in Eq. EJ) corresponds to the motion of the junction in the absence 
of beads, with an effective junction mass related to its own mass and the bead mass, with 
a center-of-mass correction for finite quark masses. The second part (terms three and four) 
is the independent motion of the beads on the three triads with respect to a fixed junction, 
with a bead mass also corrected for center-of-mass motion for finite quark masses. There is 
also an "interaction" part where the junction interacts with the various bead modes (term 
five), where the bead modes associated with the same quark interact with each other, and 
where the modes on triads corresponding to different quarks interact with each other (term 
six). Note that these bead self-interactions (term six) vanish for infinite quark masses. 

Because the quarks move with fixed relative positions only to maintain the center of mass 
position in the presence of a moving junction and beads, there are no quark kinetic terms 
in this string Hamiltonian, and there is no sense in which the quarks acquire mass from the 
beads, i.e. constituent quark masses are not derived from current quark masses. 

Note that the model predicts that for bead motion in the small-oscillations approximation 
the potential has the customary, and phenomenologically important sj], linear potential 
term &^i=i h ( see Eq. E|). For large ij, where the small-oscillations approximation becomes 
exact, the potential is just the linear term expected in any string model. This potential is 
a prediction of the model, not an ansatz. In the numerical work that follows, the small- 
oscillations approximation for the junction will be removed to yield the potential energy in 
the absence of beads, i.e. 



/ 



A/3 



h + x 



h 



+ y +z z 



(^-x) +(|+^) 2 + * 2 + V* 2 + (*3-2/) 2 + 2 2 



(13) 
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B. Hamiltonian for 9j > 120° case 



The correct expression in the small-oscillations approximation for the flux Hamiltonian 

for the case 0j > 120° can be obtained (see Appendix |BJ) by setting l 3 = N 3 = in -£% ux , 

M e g, and m e g in Eqs. EJ El and HI and g| m , q\ m = in -£% ux in Eq. El and changing the 

junction potential energy in Eq. [7| to 

2 n 

yjm.osc. =bJ2h + 2% cos — + by/x 2 + y 2 + z 2 
i=l 



6 

+ 2 



2 2 @J , 2 ■ 2 ®J , 2 

x cos — + ?/ sm — + z 



i 1 1 0j r i i , , 

/; + ^f- 2 ^ sin y cos Yi^-^^- (14) 



Note that the third term above, which is the length of the third triad when the junction has 
moved, cannot be expanded in the small-oscillations approximation. The vectors e l become 

e = I -cosy,smy,0 1 , e = I - cos —, -sm —, J . (15) 

In the numerical work that follows, the small-oscillations approximation for the junction will 
be removed to yield the potential energy in the absence of beads, 



Vj = b I W (h sin % + x) 2 + (k cos % + y) 2 + z 2 



+ x l(l 2 sm^--x) 2 + {l 2 cos^- + y) 2 + z 2 + ^x 2 + y 2 + z 2 ) . (16) 



C. Approximate flux Hamiltonian 

It will now be demonstrated that the interaction terms in Eq. EJ (terms five and six) 
give a minor contribution in the small-oscillations approximation. The free parameters in 
the model (and the values initially used for the numerical simulation) are the string tension 
(0.18 GeV 2 ), the ratio of the junction and bead masses (1), and the quark masses (0.33 
GeV). The simulation is performed with one bead between each quark and the junction, 
and the quarks at first form an equilateral triangle with the lengths of the triads given a 
typical value of 2.5 GeV -1 . 

For the purposes of this demonstration, the problem is first solved numerically without 
approximations, by solving the classical Euler-Lagrange equations of motion rather than 
using quantum mechanics. This solution should provide a good indicator of how the mode 
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frequences with and without the interaction Hamiltonian compare. The mode frequencies 
parallel and perpendicular to the plane are (in GeV) 

Parallel 0.607 0.607 0.924 1.08 1.08 

Perpendicular 0.828 0.924 0.924 1.37 

where the lowest frequencies have been identified in bold face. If the interaction Hamiltonian 
is set to zero, the number of modes does not change, as the number of degrees of freedom 
(three junction coordinates and two transverse bead coordinates per bead) is unchanged. 
Since the junction and bead degrees of freedom become uncoupled, it is possible to identify 
modes involving junction motion and those involving bead motion. The mode frequencies 
(in GeV) corresponding to the junction motion (bold) and bead vibrations are 

Parallel 0.614 0.614 1.00 1.00 1.00 

Perpendicular 0.869 1.00 1.00 1.00 

The similarity of the frequencies of the lowest energy modes in this approximation to those 
arising from the full Hamiltonian (a deviation of 1% in the case of the modes with motion 
parallel to the QQQ plane, and 5% for the mode perpendicular to the plane) shows that 
for these lowest energy modes the interaction Hamiltonian can be safely neglected. In 
retrospect, the reason for this is because of the choice of physically appropriate coordinates 
for the problem, i.e. the junction coordinates and the coordinates of the beads transverse 
to the triads joining the junction to the quarks (see Fig. |7j). 

To ensure that this result is not dependent on this choice of QQQ configuration or the 
parameters in the Hamiltonian, the parameters were varied independently around the central 
values used above. Quark masses up to the charm quark mass of 1.5 GeV were used, the 
ratio of the junction to the bead mass was taken up to 10, and the triads were given 
lengths from 0.5 — 5 GeV -1 , and cases with unequal lengths were tested. The percentage 
difference between all nine mode frequencies calculated with the full Hamiltonian and with 
the interaction terms neglected is shown for selected parameters in Table |H] The largest 
error for the two lightest parallel modes and the lightest perpendicular mode, shown in bold 
face, is 7%. This demonstrates that, to a good approximation, the dynamics of the three 
lowest frequency modes can be simplified to junction and bead motion which are independent 
of one another. The bead motion on various triads, and bead motion in various modes on 
the same triad are to a good approximation independent of each other. 
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TABLE II: Percentage differences between all nine mode frequencies obtained in the one bead per 
triad problem for the full Hamiltonian (Eq. [SJ) and the Hamiltonian with the interaction terms 
(the last two terms in Eq. |5J) neglected, for selected parameters. The first group of five percentage 
differences is for parallel mode frequencies, the last four for perpendicular mode frequencies. Within 
each group the percentage differences are displayed from left to right for in ascending order of mode 
frequency. Quark masses are in GeV and triad lengths in GeV -1 . 



Mi 


M 2 


M 3 


mj 
m h 


h 


h 


h 


Percenta 


5e difference (%) 


0.33 


0.33 


0.33 


1 


2.5 


2.5 


2.5 


1 


1 


8 


8 


8 


5 


8 


8 


37 


1.5 


0.33 


0.33 


1 


2.5 


2.5 


2.5 


2 


2 


5 


11 


11 


7 


5 


5 


40 


0.33 


0.33 


0.33 


10 


2.5 


2.5 


2.5 








2 


1 


1 





2 


2 


5 


0.33 


0.33 


0.33 


1 


0.5 


2.5 


2.5 


1 


6 


9 


2 


4 


7 





10 


11 


0.33 


0.33 


0.33 


1 


5 


2.5 


2.5 


1 


1 





2 


7 


3 


2 


8 


28 


0.33 


0.33 


0.33 


1 


0.5 


2.5 


5 


1 


1 


2 
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The frequencies can be followed from the non-interacting interactions are turned 

on, and level crossing does not occur. Hence mode frequencies for the fully-interacting 
Hamiltonian can be uniquely associated with modes frequencies obtained with interactions 
neglected. The lowest frequency is always associated with the lowest junction excitation. 
However, the next lowest frequency can be associated with the second junction excitation 
or with a bead excitation along a triad when the QQQ configuration is asymmetric. This 
work focuses on the lowest lying excitation of the flux configuration, always corresponding 
to junction motion, but it should be kept in mind that the next lowest hybrid baryon may 
involve bead excitation. 

The equal-mass three-bead model is unrealistic when one of the triads is short, since the 
mass of the bead on the short triad is not representative of the energy stored in the triad on 
which it lies. An alternative model has been considered where the bead mass is taken to be 
proportional to the length of the triad, but the sum of the bead masses and junction mass 
still equals the energy stored in the Y-shaped string configuration (Eq. El with Ni — 1). In 
this model it is found that the low-lying frequencies for the full Hamiltonian are very similar 
to the former model, with similar small errors induced by neglecting the interaction terms 
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in the Hamiltonian of Eq. 



D. Analytical solution of the flux Hamiltonian 

For the Hamiltonian in Eq. |5] (the 6j = 120° case) the last two terms have been shown 
to be negligible for the lowest frequency in the small-oscillations approximation, in the case 
where there is one bead on each triad, and these terms are neglected in what follows. What 
follows is, therefore, based on the approximate 6j = 120° Hamiltonian 

1 1 3 N 4- 1 

#flux = 2 M efF 2 + Vj + ^effE^f- E M\m? + 

i=l m=l 

i=l m=l 

with Vj given by Eq. [T31 Note that since the ql, q]_ m are defined with respect to the line 
from the junction to the quarks, they depend implicitly on r. The corresponding Hamiltonian 
for 9 j > 120° is obtained by setting l 3 = N 3 = in M e g and m e g in Eqs. El and |H1 and 
restricting the summation in Eq. El to i = 1,2, with Vj given by Eq. 

Since Eq. [T7| is diagonal in the coordinates of the beads, the last two terms in Eq. H7t 
corresponding to the bead energies, can be replaced by their ground-state harmonic oscillator 
energies 

1 3 / hi- Ni 

^- 2 M^ + V 3 ^^^E<, (18) 

which is summed over the two polarizations possible for each bead vibration. The sum of 
the frequencies is 



m=l H blU 4(Wi+l) 

If the number of beads is taken to infinity, the dependence of the Hamiltonian on the 
(unphysical) number of beads can be removed. The part of the Hamiltonian in Eq. fTHl 
arising from the beads becomes infinite when Ni — > oo. To avoid infinite energies, the bead 
separation regularization parameter (analogous to the lattice spacing) a = k/(Ni + 1) is 
fixed when iVj — > oo. This is consistent with the flux-tube model philosophy that a cannot 
be chosen arbitrarily small, since that would lead to the breakdown of the strong coupling 
expansion of the Hamiltonian formulation of QCD, from which the model is motivated [13]. 
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Taking Ni — > oo, while keeping a fixed, the Hamiltonian becomes 



3 



5 1 7T 



i=l 



7ra 2 a 12/, 



+ 0{a 2 ) + ■ 



(20) 



with 



8=1 



where is now independent of the (unknown) ratio mj/rrib. Since N{ — > oo with a 

fixed, the part of the Hamiltonian in Eq. ED arising from the beads is valid only in the limit 
Li — > oo. 

The part of the Hamiltonian in Eq. ED arising from the beads contains a linear term 
Ali/ira 2 and a constant term — 1/a, which are regularization-scheme dependent terms in 
the "self energy" of the string system. As explained in Ref. Q|, the linear term should 
be regarded as a contribution that renormalizes the bare string tension b to its physical 

r~i 

value. The constant term is three times larger than the constant term found for mesons . 



The Luscher |2J] term — 7r/12/j is regularization-scheme independent and finite, and can 
be regarded as a prediction of the flux-tube model, although it is insignificant at large Zj, 
where its derivation is valid. The Luscher term arises in relativistic string theories |2J| in 
the limit k — > oo. In this limit the string excitations in our model coincide with relativistic 
string theories. There exists strong lattice-QCD evidence, through the study of "torelons" , 
of the existence of the Luscher term in QCD 25]. The form of the Luscher term should be 
contrasted with that of a Coulomb term, with the former depending inversely on the triad 
lengths, and the latter inversely on the distance between two quarks. 

The Hamiltonian for the 8j > 120° case is obtained by setting I3 = in in Eq. 
and restricting the summation in Eq. [201 to i = 1, 2 with Vj from Eq. 

If the redefined adiabatic approximation was not made, i.e. if the calculation was not 
performed in the center of mass frame of the entire system with the distances between 
the quarks fixed, then = &^Zj/3. The correction from center-of-mass motion in 

Eq. |^ substantially reduces the effective mass of the junction. It is shown below that the 
excitation energies of the junction are proportional to (i\/p^)~5 (see Eqs. IHDTl31|) . and so 
with typical QQQ configurations the junction excitation energies are 1 to 2 times larger in 
the redefined adiabatic approximation than in the adiabatic approximation. This underlines 
the importance of working in the CM frame. 
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The ground state bead configuration that solves the Hamiltonian of Eq. El is 

3 Ni r— 

8=1 771=1 

ex P " 2 E E ^eSbkm + 1) ^ [(gf,J 2 + fej 2 ] , (22) 

V 1=1 771=1 J 

in the limit N% — > oo (with a fixed) that is used to express the Hamiltonian in Eq. I2H The 
corresponding wave function for 9j > 120° is obtained by restricting the products in Eq. 1221 
to i = 1, 2, and setting Z 3 = in m e g in Eq. |H1 



E. Analytic small-oscillations solution to the junction Hamiltonian for 6j = 120° 
Define 

r i ~~ r 2 x r i + r 2 — 2r 3 p • A 

The six Jacobi variables p, A consist of (1) four variables which specify the positions of the 
quarks in the QQQ plane, p = \p\, X = |A|, and 9 P \ (the angle between p and A) and 4> p 
(the angle between p and the space-fixed x-direction) , and (2) two polar angles 9, and <p 
which specify the orientation of the vector p x A which lies perpendicular to the plane. The 
variables p, A and 9 P \ are rotational scalars. They are related to the triad lengths l\, I2, and 
Is by the relations 



1%—N (p 2 + pX sin 9 pX + V3pX cos 9 P \) l 2 = M (p 2 + pX sin 9 P \ — V3pX cos 



P 2 



h—N (-y + ^A 2 + P A sin pA ) A/" -1 = a/ ~(p 2 + A 2 ) + 3pA sin 9 pX . (24) 

In the remainder of this section junction motion is treated in the small-oscillations ap- 
proximation for the 9 j = 120° case. The junction Hamiltonian is |M^r 2 + Vj m - osc -, where 
the potential is from Eq. [7| Junction motion in (x, y) and out (z) of the QQQ plane are not 
coupled in Vj m-OSC- , so motion along the z direction is one of the vibrational modes of the 
junction. One way to define the z direction in terms of the positions of the quarks is by the 
normalized vector 

P x ^ ^ , 

= "-fesr- (25) 

where a z denotes a sign that will be specified later. Note that z motion of the junction 
motion occurs along the direction of the vector p x A, but there is no physical reason to 
prefer one sign a z over the other. 
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The in-plane part of the Hamiltonian can be diagonalized in terms of the normalized 
eigenvectors j^] 



N(h,l 2 ,h) 



1 , 1 



h 2 y i l + T 2 )±v/ ^ ll,k,h) 



^VSA 1 



where 



and 



1 1 



s(/i, h, h) = To + 



il ^2 ^3 ^1^2 ^1^3 



/■/ 



>0, 



2 '3 



N(l u k, h) = \/2js(h, l 2 , h) ± (I - 1(1 + 

V '3 ^ *1 '2 



3j- 



(26) 



(27) 



(28) 



and a sign a± is included for the same reasons as above. The vectors fj z , fj + , and fj _ can 
be verified to be orthonormal vectors. The junction Hamiltonian can now be written as 

3 



M°° f 2 -I- T/ S1 



i=l 



; M eff [($+ ■ + (V- ■ + (Vz ■ r) 2 + ^+0?+ • r) 2 + u 2 _ (fj_ ■ r) 



uj 2 z {f, z ■ if] ,(29) 



where the vibrational frequencies are given by 

b , 1 



and 



2M eff 



1 u 

M^ y h + l 2 + lj 



111 \ 
T + - r + - r ±V s (h,h,h) ■ 
n h h J 



(30) 



(31) 



Note that the out-of-plane mode is always more energetic than the in-plane modes, since 
uo z > u)±, and that the in-plane modes have uo_ < u + , with degeneracy only when l\ = l 2 = I3. 

Solving the Schrodinger equation corresponding to Eqs.[2IIl and I2H1 yields the ground-state 
energy, corresponding to the adiabatic potential for the quark motion in a conventional 
baryon, of 

3 3 



4L 



1 



7T 



7TGr 



a 12L 



0(a 2 ) + 



(32) 



V B (h,k,h) =bJ2h + ^(LU + +LU_+LJ z ) 

1=1 i=l L ' 

Junction excitations in the fj_, r) + , or f) z directions yield adiabatic potentials for the quark 
motion in different low- lying hybrid baryons, denoted H\, H 2 , and H3, ordered from least to 
most energetic. The hybrid baryon string energy (adiabatic potential) is that of the baryon 
in Eq. EH with the term c<j_, uj + or uj z added for Hi, H 2 or H 3 hybrid baryons, respectively. 
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Note that these results neglect the junction-bead and bead-bead interactions, which has 
been demonstrated to be a good approximation only for the lowest energy (a>_) mode. 

It is intriguing to note that the baryon potential in Eq. E21 serves as an analytical form 
to which the infinitely-heavy quark potentials calculated in lattice QCD can be fitted as 
a function of Zj. Furthermore, the ij dependence of the various hybrid baryon potentials 
are predicted. Potentials were also predicted in Ref. j^J. Comparisons to lattice results 
would be instructive at large U, for which Eq. l3*2l was derived. The physical string tension 
is b — 4/ (ira 2 ). The constant — 1/a term is regularization-scheme dependent, and hence not 
physical. Indeed, lattice calculations [15| find the constant term regularization dependent, 
and proportional to 1/a. The remaining terms do not depend on either b or a when Mj — > oo, 
which is the limit in which lattice QCD potentials are evaluated, noting that b/M^ is 
independent of b in this limit. 

The normalized flux wave function of the baryon is 



7T4 



roc 



tt B (r) = — ^ 3 exp ^ ^(^+(r7 + • r) 2 + u„(r)_ ■ if + u z {t) z ■ r) 2 ) \ $, 



(33) 

where the bead wave function $ in Eq. |221has been incorporated. For the Hi, H 2 , and H 3 
hybrid baryons, respectively, the normalized flux wave functions are ^/^(r) 



2M^< 



f i _ ■ r * B (r), VSMrfK ^+ " r ^s(r), J2M^j g fj z ■ r * B (r). (34) 



F. Junction motion away from the small-oscillations limit 

Eq. ^1 cannot be expanded in the small-oscillations approximation to the junction mo- 
tion. Without this approximation, i.e. where Vj in H^ in Eq. 1201 is taken from Eqs. 
[TBI and [TBI for the cases 6j = 120° and 9j > 120° respectively, the eigenfrequencies and 
eigenvectors cannot be solved for analytically. 

The variational principle is used to separately minimize the expectation value of the 
Hamiltonian H^ by solving the Schrodinger equation for the conventional baryons and 
hybrids Hi using the ansatz simple-harmonic-oscillator wave functions in Eqs. E21 and 13*41 
The calculated energies are upper bounds for the true energies, according to the Hyleraas- 
Undheim theorem. The parameters of the ansatz wave functions no longer have the values 
that they had in the small-oscillations approximation, but need to be fitted. For example, 
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the directions i) z , fj_, and r) + are no longer given by Eqs. E^andES but will be fixed by 
the variational principle. 

Note that (Eq. |2*U|) is even under the discrete transformation z — > —z since Vj in 

Eqs. El and El only depends on z 2 . This implies that the wave functions should be either 
odd or even under z — > —z. But since the wave functions are assumed to be of the form in 
Eqs. Hn3and|jMJ it is not difficult to show that this implies that one of the junction vibrational 
modes, corresponding to rj z , is always perpendicular to the QQQ plane. Note that fj_, fj + , 
and fj z are required to be orthonormal in order to obtain orthonormal hybrid baryon wave 
functions in Eq. El This gives four variational parameters that specify the ansatz wave 
functions: M^u_, M^u> + , M t ^uj zi and an angle that describes the ray in which fj_ lies in 
the plane relative to the x-direction defined in Fig. The minimization is carried out with 
respect to these four variables. 

III. FLUX SYMMETRY 

A. Quark label exchange symmetry 

Denote by P 12 , P13, and P23 the permutations which exchange the labels of the quarks. 
Except for color, quark-spin and flavor labels, which will only be of interest later, exchange 
symmetry affects only position labels. Under such quark label permutations the positions of 
the quarks are exchanged, e.g. P\i exchanges r x *^ r 2 , but note that variables that are not 
functions of the are unaffected. Since the physics does not depend on the quark position 
labelling convention, the flux Hamiltonian given by Eq. |S] should be exchange symmetric. As 
the equilibrium junction position r e q in Eq. ^is invariant under the P^-, and Zj = rj — r e q, 
it follows that Zj «-> Zj under Pjj. Also, since the number of beads on the i-th triad is 
Ni = U/a — 1, then Ni <-> Nj under P^-. The potential Vj in Eq. fTHl can be written in the 
manifestly exchange-symmetric form 

3 

V J = b J £\n-r\, (35) 

i=l 

noting that the junction position r is not determined by the positions of the quarks. This 
establishes that all quantities in the flux Hamiltonian in Eq. 03 for the 8j = 120° case are 
invariant under exchange symmetry transformations. 
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Since the flux Hamiltonian is invariant under exchange symmetry, it is clear that energies 
(or adiabatic potentials) that are solutions of the flux Schrodinger equation are also exchange 
symmetric. This is explicit for the frequencies in Eqs. ISHandHni and the potential in Eq.l32l 

By the same arguments as above it can be shown that in the 8j > 120° case, the Hamil- 
tonian (Eq. I14j) is invariant under P 12 . 

Since the Hamiltonian is exchange symmetric, the commutation relations [Hq ux , Pij] = 
hold. This implies that the wave functions of conventional and hybrid baryons have 
to represent the permutation group S3. Possible representations are the one-dimensional 
symmetric and antisymmetric representations, and the two-dimensional mixed-symmetry 
representation. Since the baryon and each of the hybrid baryons Hi have different [s^.] flux 
energies V(li, l 2 , h), where H^ ux ^/ = V(li, l 2 , £3)^, each of the four wavefunctions ^/_e(r) 
and ^^(r) have to belong to a one dimensional representation, as they cannot mix with 
each other under any of the permutations Py. This implies that \& #(r) and \I/#.(r) are either 
totally symmetric or antisymmetric under quark label exchange. 

In the baryon wave function of Eq. ESI the quantities 00 u z , and are exchange 
symmetric, so that the factor before the exponential is invariant. The bead wave function $, 
given in Eq.|221 is also invariant. Since ^^(r) is either exchange symmetric or antisymmetric, 
the exponential function in the junction coordinates must be either exchange symmetric 
or antisymmetric. The second possibility is untenable since the exponential function is 
always positive. Hence, the baryon wave function ^^(r) is totally symmetric under exchange 
symmetry. 

Consider the hybrid baryon wavefunctions in Eq. IS31 The above implies that f)_, r) + , and 
r] z are either totally symmetric or totally antisymmetric under exchange symmetry, since r 
is independent of the quark labels. It is shown in Appendix that both possibilities are 
explicitly realizable. This implies that for each of the hybrid baryons Hi, there is a degenerate 
pair of totally symmetric (S) and totally antisymmetric (A) wave functions, denoted by Hf 
and Hf. 

The preceding argument assumed that \1/ has the form in Eqs. ESI arid H2J which applies 
only in the 9j = 120° case with small junction oscillations. However, as was discussed 
in section Hi F| in the more general case where small junction oscillations are not assumed, 
ansatz wave functions of the form in Eqs.|SSl an d|Mlare used, so that the preceding arguments 
regarding exchange symmetry remain valid. 
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The above arguments can be repeated to show that in the 6j > 120° case, the baryon 
ansatz wave function \l/s(r) is invariant under Pu, and that the hybrid baryon ansatz wave 
functions ^^(r) and fj_, fj + , and fj z are either odd or even under P 12 . 

B. Parity 

The operation of the inversion of all coordinates, or parity, is a symmetry of the flux 
Hamiltonian. It follows that i) z in Eq. 123 is even under parity, since p — > — p and A — > —A 
under parity. If a z — 1 this follows trivially, and if a z is given by Eq. IC21 it follows because 
Eq. IC2I is invariant under parity. 

The U remain invariant under parity since they are lengths, but the Z; are odd under 
parity. From the definition of the r) ± in Eq. EH1 and the definition of x and y in terms of 
the li in Eq. it follows that the fj± are odd under parity. The sign a± is invariant under 
parity (see Appendix EJ) • This argument is so far valid only when fj± is given by Eq. EEl 
applicable for the 6j = 120° case in the small-oscillations approximation. However, for the 
ansatz variational wave functions in section Hi Fl the f)± lie in the QQQ plane and so must be 
linear combinations of p and A with coefficients which are functions of the parity-invariant 
variables p, A and 9 pX , so the r) ± remain odd under parity. 

Since the position r of the junction is a vector, it is odd under parity. It follows that the 
baryon wave function \I/,b(i") in Eq. ESI is invariant under parity. The hybrid baryon wave 
functions in the QQQ plane, i.e. ^/^(r) and \I/# 2 (r) in Eq. are even under parity, while 
\1/ H :j { r ) is °dd under parity. These results also obtain for 6j > 120°. 

In summary, flux wave functions of baryons and H\ % % hybrid baryons are even under 
parity, while the H 3 hybrid baryon flux wave functions are odd under parity. 

C. Chirality 

Reflection in the QQQ plane, or "chirality" j^l, is generally a symmetry of the flux wave 
function in the adiabatic approximation, since the physics does not distinguish between 
above and below the QQQ plane. The relevant group consists of the identity and reflection 
transformations. In this approximation the flux wave function can be classified according 
to its eigenvalue under reflections in the plane spanned by the three quarks, which is the 
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chirality ±1. 

In the flux-tube model this reflection takes z — > —z and q]_ m — > —q± m . The most general 
Hamiltonian derived in this work, Eq.^1 is invariant under this reflection transformation, as it 
must be. The baryon and hybrid baryon wave functions in Eqs.|2Sland|^]are eigenfunctions 
of the reflection transformation. The baryon and "planar" hybrids (#1,2) have chirality 1, 
and the "non-planar" hybrid H3 has chirality -1. Hence the chirality formally allows us to 
clearly distinguish "planar" and "non-planar" hybrids, even for more general solutions of 
the Hamiltonian than those in Eqns. ESI and 

In adiabatic lattice QCD, exchange symmetry, parity and chirality should classify the 
(hybrid) baryon flux wave functions. These properties are sometimes called the "quantum 
numbers of the adiabatic surface" . 



IV. QUANTUM NUMBERS 

A. Orbital Angular Momentum 

For every set of quark positions the potential in which the junction moves is anisotropic, 
which means that the solutions of the Schrodinger equation for the junction motion do not 
have definite orbital angular momentum or its projection. However, in the absence of the 
adiabatic approximation the combined wavefunction of the quark and junction motions must 
be a state of good angular momentum. 

It is possible to determine the angular momentum character of the variational wavefunc- 
tions which minimize the flux energy for a given set of quark positions. The probability 
of overlap between an isotropic S*-wave harmonic-oscillator state with frequency u and the 
baryon flux wavefunction ^s(r) of Eq. IH31is 

m = 0) B |(>Mr)|000)| 2 = V^-- ^. (36) 

(cj_ + Uj) [UJ + + UJ) (UJ Z + UJ) 

and that of an isotropic P-wave harmonic-oscillator state with the flux wavefunction of the 
lightest hybrid (r) of Eq. EH is 



Px(i = l) S £ |(tt Hl (r)|01M)| a = ^if^^ ~ y (37) 
M ±-f +1 (w_ + uj) 3 (uj + + w)(w z + UJ) 

Once the energies of the ground and first excited states of the flux have been independently 
minimized in the variational calculation described below in Sec. El the calculated values of 
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TABLE III: Flux energies and angular momentum probabilities calculated using Eqs. EHl and 1371 
with y/ujM^ = 0.4 GeV, for four quark configurations. Here £V M t = 0.99 GeV, b = 0.18 GeV 2 , 
the triad lengths l{ are given in GeV"" 1 , and energies are in GeV. Here Eq and Pq(1 = 0) are the 
energy and S- wave probability of the flux ground state, E\ and Pi(l = 1) are the energy and 
P-wave probability of the first excited state of the flux. 



h 


h 


k 


Eq 


Foil = 0) 


Ex 


Pi(l = l) 


2.5 


2.5 


2.5 


1.09 


0.995 


1.76 


0.997 


2.5 


2.5 


0.5 


1.42 


0.999 


2.18 


0.998 


2.5 


5.0 


0.5 


1.18 


0.993 


1.80 


0.998 


0.5 


0.5 


10.0 


1.30 


0.986 


1.92 


0.998 



co>_, oj + , and u z can be used to find these probabilities. The result of these numerical studies 
is shown for sample quark configurations in Table IIII1 It is clear that the ground state of 
the flux is in an almost exclusively angular momentum zero state, so the orbital angular 
momentum of the baryon is that of the quark motion. 

Table IIHI shows that variational calculations result in flux wave functions ( r ) which 
are to better than 99% a linear combination of Yu(r) and Yi_i(r). An alternative argument 
is given here that the angular momentum of the flux in the lowest-lying hybrid baryons (Hi) 
is predominantly unity. The flux wave function in Eq. EH of the lightest (Hi) hybrid baryon 
is proportional to rj_ ■ r, where i)_ lies in the plane of the quarks. If the exponential in 
Eq. |^]was spherically symmetric, it would be strictly true that ^^(r) was proportional to 
a linear combination of VJi(r) and Y/_i(r), with the junction position r denned relative to 
a (body) z-axis perpendicular to the quark plane. 

Further numerical studies, described below, have shown that the least energetic motion of 
the quarks in the Hi adiabatic potential has the quark angular momentum L q = 0, the next 
highest L q = 1, etc. Furthermore, there is a substantial cost in energy to increase the quark 
angular momentum in the Hi hybrid potential, so the total orbital angular momentum of 
the lightest hybrid baryon is unity. 

In order for the combined flux and quark orbital angular momentum to have a definite 
value (unity), in principle the components with orbital angular momentum other than unity 
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in the flux wavefunction must be combined with quark motion with L q > 1 to make the 
total orbital angular momentum unity. Given the negligible size of these components, a very 
good approximation to the energy can be found by assuming that the flux orbital angular 
momentum is unity 



B. Color 



It is important to note that the wave function of the (hybrid) baryon has both a color 
and a flux sector, which are separable. This is because color is a separable degree of freedom 
in quantum chromodynamics, which labels the quarks and flux lines. This is described by 
the color sector of the theory. The flux sector, on the other hand, concerns the dynamics of 
the flux. In the bag model and large N c limit the same separation occurs, where the octet 
color of the gluon is combined with that of the quarks, and the spatial motion of the gluon 
is treated separately 0,0]. 

In the flux-tube model the color structure of a hybrid baryon is motivated by the strong 
coupling limit of the Hamiltonian formulation of lattice QCD [28]. Here, the quarks are 
sources of triplet color, which flows along the triad connecting the quarks to the junction, 
where an e-tensor neutralizes the color. The color wave function is hence totally antisym- 
metric under exchange of quarks for both the conventional and hybrid baryon. In the bag 
model and in the large N c limit the color structure of a hybrid baryon is very dif- 
ferent. This color structure is critical for the correct exchange symmetry properties of the 
conventional and hybrid baryons, and hence the structure of the wave function. 



C. (Hybrid) Baryon Wave Functions 

The energy of the quarks in the potential given by the flux energies is found by expanding 
the quark wave function in a basis with well defined orbital angular momentum L q and 
projection M q , made up from orbital angular momenta l p and l x in the coordinates p and A 
respectively, 

(p,X\n p l p n x l x ;L q m q ) = N nplp R np i p {p) Af nx i x R nx i x (\) 

x C(l p m p l x m x -,L g M g )Y lpmp (n p )Y lxmx (n x ), (38) 

m p m x 
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where M n i is a normalization factor, and the Clebsch- Gordon coefficient combines spherical 
harmonics with orbital angular momentum l p and l\ to form a state with orbital angular 
momentum L q . Here the R n i are orthonormal and complete functions in the radial co- 
ordinate, where n = 0,1,2... denotes the radial quantum number, which are taken to be 
three-dimensional harmonic oscillator radial wave functions, i.e. Laguerre polynomials. It is 
easy to show that the wave functions in Eq. EH1 form an orthonormal (in all six labels) and 
complete set. In Eq. EB]a formal notation is used where the wave function is defined as the 
overlap of a state \n p l p n\l\; L q M q ), characterized by the quantum numbers indicated, with 
a position state |p, A). 

A linear combination of the states in Eq. EH can be used to form a general eigenstate 
of quark orbital angular momentum L q and projection M q , denoted by \nL q M q ), where n 
denotes the radial quantum number. The corresponding wave function is 

(p,X\nL q M q ) = c tUiJP,Mn P l P nxh;L q M q (p,\)). (39) 

The coefficients in this linear combination, and the corresponding hybrid baryon energies, are 
found by diagonalizing the three-quark Hamiltonian in the basis of Eq. EElwith the potential 
energy given by the flux energy. Note that the orbital angular momentum and spin of the 
quarks are good quantum numbers as the inter-quark potential is a spatial and quark-spin 
scalar, even in the presence of the Coulomb and hyperfine (spin-spin) interactions. 

It has been checked numerically for the adiabatic potentials found here that the lowest 
energy solutions of the Schrodinger equation for both conventional and hybrid baryons have 
L q = quark wave functions. In order to determine the color, flavor, quark-spin, parity, 
exchange symmetry and chirality quantum numbers of these states, it is sufficient to consider 
the (p, A|0000; 00) component of the L q = wave function in Eq. H2B as these quantum 
numbers must be the same for all components of the wave function. 

From Eq. HH (p, A|0000; 00) equals 

a 3 f a 2 1 
A/- 2 i?oo(p) ^oo(A) Y 00 (Q p ) Y 00 (Q X ) = -3-exp j-y(p 2 + A 2 ) j , (40) 

where a is a parameter that characterizes the Laguerre polynomials. This is obviously even 
under parity and is totally symmetric under exchange since p 2 + A 2 is exchange symmetric. 
Since the parity is unaffected by the color, flavor, and quark-spin wave functions which will 
multiply this L q = spatial wave function, the parity is determined by that of the flux 
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TABLE IV: Quantum numbers of ground state L q = flux-tube model (hybrid) baryons for the 
lowest flux-tube surfaces B (the conventional baryon) and H\ (the lightest planar hybrid baryon). 
In the absense of spin-dependent forces all ground states corresponding to a given flux-tube surface 
[both symmetric (S) and antisymmetric (A)] are degenerate. Here L is the total orbital angular 
momentum of the quarks and the flux, N or A denotes the flavor, S is the spin of the three quarks, 
J = L + S is the total angular momentum, and P is the parity. Low-lying hybrid baryons in the 
bag model contructed with a transverse electric gluon (corresponding to the surfaces Hf and H±) 
are also shown Ifjl. 
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wave function given in section UlI Bl The parities of the low-lying hybrids are displayed in 
Table ITVl 

The quark-spin (x) x flavor (<p) wave function can be made totally symmetric for quark- 
spin | x flavor A, using the product of symmetric factors xf/2^A> an d f° r quark-spin | x fla- 
vor A/, using the linear combination 2^ of mixed-symmetry factors (x 1 /20Jv P +Xi/2^iV A )/v / 2- 
It can also be made totally antisymmetric for quark-spin | x flavor using the linear com- 
bination of mixed-symmetry factors (Xij^N* ~ X^'Pn'')/^- 

Since quarks are fermions, the combined color, space, quark-spin, flavor and flux wave 
function should be totally antisymmetric under exchange symmetry. Since for L q = 
baryons and hybrid baryons the color and space parts are totally antisymmetric and sym- 
metric, respectively, the flavor x quark-spin x flux part must be totally symmetric. 

For baryons the flux wave function is totally symmetric with orbital angular momentum 
zero, and so their quantum numbers are exactly as they were in the conventional quark 
model with an assumed static confining potential between the quarks. As an example, 
the quantum numbers of the non-strange L q = ground states are shown in Table IIV1 
The L q = hybrid baryons Hf have a totally symmetric flux wave function, and so the 
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quark-spin x flavor structure is the same as for the corresponding L q = baryons, i.e. 
the symmetric products above with quark-spin 1/2 for nucleons, and quark-spin 3/2 for A 
states. For L q = hybrid baryons with a totally antisymmetric flux wave function, Hf, the 
quark-spin x flavor wave function must be totally antisymmetric, and the only possibility 
is the antisymmetric product with nucleon flavor and quark-spin 1/2, as shown in Table HV1 

Chirality is a reflection in the QQQ plane, and hence only affects the flux part of the 
wave function, so that its values are those given in section UlI CI 

For the lightest L q = (hybrid) baryons the total orbital angular momentum is that of 
the flux. This gives L = for low-lying conventional baryons, so that J = S. Since L = 1 
for the low-lying Hi hybrid baryons, J = | or | for S = ~, and J = \ ,|, or | for S = |, as 
shown in Table HVl 



V. HAMILTONIAN FOR THE QUARK MOTION 

A phenomenological form is used here for the quark Hamiltonian which is fit to conven- 

n 

tional baryon spectroscopy in Ref . |8| . In the case of hybrid baryons the difference between 
the adiabatic potential found from numerical calculation of the energy of the ground state 
and the first excited state of the flux is added to the quark Hamiltonian. 
The quark Hamiltonian has the form 

3 , 

Rqqq = E v p * 2 + M ? + E v S oul + E ^ cont + v(h, h, h), (41) 

i=l i<j i<j 

where P« is the momentum operator of the i-th quark, Mj = 0.22 GeV for light quarks, 
b = 0.18 GeV 2 , and the Coulomb potential V^ oul and hyperfine contact potential V? ont have 
the same form as in Ref. [sj]. The justification for adopting this form of the Coulomb and 
hyperfine contact interaction is outlined in Sec. IV Bl below. For the conventional baryon the 
adiabatic potential Vb(Ii, h, h) also has the form b U adopted in Ref. In what follows 
the numerical calculation of the form of the adiabatic potential V for the lightest hybrid 
baryon is outlined. 
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FIG. 8: Difference Vh 1 — Vb of the hybrid and conventional baryon adiabatic potentials for p = 6.2 
GeV -1 , as a function of A and 6 p \. 

A. Numerical adiabatic potentials 

As it is not possible in all quark configurations to derive the adiabatic potential of baryons 
and hybrid baryons in the small-oscillations approximation to the junction motion given by 
Eq. E21 a numerical calculation is used to find the flux energy, which is part of the potential 
energy for the quark motion, for all quark configurations. As discussed previously, the linear 
term in Eq. EH which arises from the bead motion is regularization-scheme dependent, and 
can be absorbed into the physical linear term in the potential. Also, there will be no need to 
consider constant and Luscher terms in this section as they are identical for the conventional 
and hybrid baryons. 

The procedure of numerically evaluating the hybrid baryon potential is as follows. For a 
large set of quark configurations {I1J2J3}, the Schrodinger equation 



is solved variationally for the wave function ^s(r) and energy Vb(Ii, h, h) of the ground state 
of the junction Hamiltonian, as described in section Hi Fl using the ansatz wave function in 
Eq.Eni The lowest-lying hybrid (Hi) baryon potential is found by solving Eq. l42lvariationally 
for the wave function (r) and energy Vh 1 (h, h, h) of the first excited state of the junction 
Hamiltonian, using the first ansatz wave function in Eq. 0U This corresponds to junction 
motion in the QQQ plane, and is used because the analytical solutions in section III El 
suggest that the lowest hybrid baryon energy can be described by such junction motion. 




(42) 
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FIG. 9: Difference Vh 1 — Vb oi the hybrid and conventional baryon adiabatic potentials for 
9p\ = vr/2 as a function of p and A. 

The minimizations for the baryon and lightest hybrid baryon potentials are carried out 
independently. It has been checked that the numerical potentials calculated here in the 
long-string limit agree with the analytic expressions derived in that limit in Sec. Ill El to 
within 2%. 

The hybrid baryon adiabatic potential is defined to be 



Selected numerical results for the difference Vh 1 — Vb of the hybrid and conventional baryon 
adiabatic potentials are plotted in Figures IBHH1 in Fig- El the potential is plotted with fixed 
p (proportional to the separation of quarks 1 and 2) and variable A (proportional to the 
separation of the center-of-mass of quarks 1 and 2, and quark 3) and 9 P \ (the angle between 
the vectors p and A), which clearly demonstrates a discontinuity in the derivative when 
the flux goes from the shape with 9j = 120° to that with 8j > 120°. In Figs. M and [TO] 
the behavior of Vh x — Vb and Vb — b J2i h when p and A are orthogonal is plotted against 
p and A. It is obvious that both the conventional and hybrid baryon adiabatic potentials 
increase when pX is small, with the hybrid adiabatic potential increasing faster. If the small- 
oscillations approximation were employed they would tend to infinity as pA — > 0. Solving 
for the energy variationally has softened this behaviour considerably. 

The value of ./ M^cu^ for the baryon is approximately in the range 0.37 — 0.5 GeV, while 



V Hl {h,l2,h) = VBihihJs) + V Hl (h,l 2 ,h) -V B (h,hj3) 




(43) 
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FIG. 10: Conventional baryon adiabatic potential without the confining potential, Vb — b^2 i h, 
for 0p\ = 7r/2, as a function of p and A. 

the hybrid baryon is ~ 0.35 — 0.48 GeV. 

B. Short-distance interactions between the quarks 

There are two important configurations of the quarks when considering the Coulomb 
interactions. These are where two quarks are near to each other and the third is distant 
(meson-like configurations), and where all three quarks are close to each other. It is possible 
to focus on the former, because the latter is atypical and contributes little to the energy of 
the baryon. 

In the flux-tube model, in the meson-like configuration a string extends from the distant 
quark to the other two quarks, i.e. the system looks like a meson with the two nearby 
quarks in color 3. The reason for this is that the long-distance picture is still appropriate 
for the distant quark, which means that its color flows along the long string, and must be 
cancelled by the two nearby quarks. The two nearby quarks are hence in color 3 for both 
the conventional and hybrid baryon. This implies that the Coulomb interactions between 
the nearby quarks are attractive, and identical in the conventional and hybrid baryon. 

This also has implications for hyperfine contact (spin-spin) interactions, the existence of 
which has recently been confirmed in lattice QCD [30j, since their form is given by the color 
representation of the interacting quark pair. The interactions are, therefore, identical for 
conventional and hybrid baryons. As seen in Table HVl above, the lightest A hybrid baryons 
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have the same quark-spin structure as the ground-state A, and both the symmetric and 
antisymmetric lightest nucleon hybrids have the same quark-spin structure as the ground- 
state nucleon. This implies that the lightest A hybrid baryons are always heavier than the 
lightest N hybrids, due to the same hyperfine interaction which makes the A heavier than 
the nucleon. 

VI. NUMERICAL MASS ESTIMATE FOR LOW-LYING HYBRID BARYONS 

The Hamiltonians in Eqs. ^2 and 1321 are evaluated using the basis of coupled three- 
dimensional harmonic oscillator wave functions in Eq. EH expanded up to at least the 
N = 7 oscillator level, where N = 2{n p + n\) + l p + l\. These matrices are subsequently 
diagonalized to yield the energies. The resulting full wave functions fEq. l3Uj) . which are linear 
combinations of the harmonic oscillator wave functions, are solutions of the Schrodinger 
equation for the quark motion in the presence of the usual Vb baryon confining potential 
and the Vh x hybrid-baryon adiabatic potential. The differences between the energies for the 
hybrid and the conventional baryon are then added to the experimental mass of the lightest 
baryon (the nucleon). 

It is interesting to determine what sets the scale of the energy difference between the 
lightest hybrid and conventional baryons in this model. This can be illustrated by examining 
the analytic solution to the junction Hamiltonian developed above, in the limit of small 
oscillations and a large number of beads. Equation EU gives the frequency uj_ of the lowest 
energy string excitation in terms of the string tension b, the effective mass of the 
junction in the limit of a large number of beads, and the lengths U of the three lines from 
the rest position of the junction to the quarks (triads). The effective mass of the junction 
is, in turn, given in terms of the same quantities and the quark masses Mj in Eq. |^ The 
scale of the energy difference is set by the same quantities in our variational calculations, 
which do not employ the small oscillations limit. 

Consistent values of the string tension and light (M u = M d ) quark mass are used in 
evaluating the flux energies, which define the adiabatic potentials in which the quarks move, 
and in the calculation of the energy of the quarks (and so hybrid and conventional baryon 
masses) in these potentials. The value of the string tension used is b = 0.18 GeV 2 , the 
typical value resulting from lattice gauge theory calculations, which is also used in quark 
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model calculations of hadron masses Mj. The light quark masses are set to 0.22 GeV, the 

n 

same value used in Ref . |8J , which also uses the relativistic kinetic energy, and Coulomb and 
hyperfine contact potentials identical to those in the current calculation. 

The triad lengths are not parameters, as in principle they are evaluated for every set 
of quark positions in order to find the potential in which the quarks move. Their average 
size does of course affect the excitation energy of the string, and since this is determined by 
solving for the motion of the quarks in the confining (and to a lesser extent short- distance) 
potential, it only depends on the quark masses and string tension. 

In the case that the hyperfine contact spin-spin term in Eq. |^is set to zero, the lightest 
L q = states have M Hl — M B = 890 MeV, giving a mass estimate of M Hl = 1085 + 890 = 
1975 MeV. Here 1085 MeV is the spin-averaged mass of the nucleon and A ground states. 
Note that this means that all of the lowest-lying Hi hybrid states in Table IIVI have this 



mass. Furthermore, the states built on this adiabatic surface with L q = 1 and L q = 2 have 
masses 2340 and 2620 MeV, respectively, showing a considerable cost in energy for orbital 
excitation of the quarks, comparable to that in the conventional baryons. Similarly, the 
lightest radial excitation built on this adiabatic surface has mass 2485 MeV, with a position 
between the L q = and L q = 2 states, as is also the case in conventional baryons. 

Including the hyperfine contact spin-spin term in Eq. |4*T1 lowers the mass of the quark-spin- 
1/2 Hi hybrid states, which coincide with the lightest N hybrids as outlined in Sec. IIV CI by 
110 MeV to 1865 MeV, and raises the mass of the quark-spin-3/2 Hi hybrid states, which 
coincide with the lightest A hybrids, by a similar amount. 

These mass predictions depend on the form of the quark Hamiltonian used here (Eq. 14 
which has been fit to the baryon spectrum. The parameters which determine the adiabatic 
potential are the string tension b and the sum of the quark masses Mj. In order to con- 
servatively estimate the error in the hybrid masses due to uncertainties in these parameters, 
a variation of Sb/b = ±10%, and S^2 i M i / ^ li M i = +50% have been considered. These 
variations change the mass prediction for the lightest hybrid by less than ±100 MeV. 

These mass estimates are substantially higher than other mass estimates in the literature, 



which are approximately 1.5 GeV in the bag mode, H and 1.5 ± 0.15 GeV in QCD sun, 
rules |4j. 

There are two crucial assumptions that were made in the early work on hybrid meson 
masses in the flux-tube model: the adiabatic motion of quarks and the small-oscillations 
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approximation for bead motion Q| . It was later shown that when the adiabatic approxima- 
tion is lifted, the masses go up, and when the small-oscillations approximation is lifted, the 
masses go down [jy] . In the present study of hybrid baryons, the adiabatic approximation 
has been partially lifted by the use of the redefined adiabatic approximation. The small- 
oscillations approximation has been fully lifted. The effect on the masses of hybrid baryons 
when the various approximations are lifted is the same as those found for hybrid mesons. 

In the numerical simulation described above, with an identical Coulomb interaction in 
conventional and hybrid baryons, the rms values a/ (p 2 ) = a/ (A 2 ) = 2.12, 2.52 GeV -1 are 



obtained for the baryon and Hi hybrid baryon sizes, respectively. The result a/ (p 2 ) = a/ (A 2 ) 
is expected since the spatial parts of the wave functions of the low-lying states are totally 
symmetric under exchange symmetry. The hybrid baryon is, therefore, 20% larger in size 
than the conventional baryon. 

It is of physical interest to obtain an estimate of the effective junction mass in Eq. |^] 
Taking the average values of a/ (p 2 ) and a/ (A 2 ) above with the quarks in an equilateral 
triangle configuration, with b = 0.18 GeV 2 , and using Eq. El the junction mass is = 
0.17, 0.20 GeV for the baryon and Hi hybrid baryon, respectively. This effective mass is 
made smaller by the center-of-mass corrections due to the quark motion in Eq. |^ This 
effective junction mass found in the flux-tube model in the refined adiabatic approximation is 
very different from the constituent gluon mass of 0.8 GeV typically employed in constituent 
gluon models Q], and partially accounts for the higher excitation energy of the hybrid in 
the present work. 



VII. DISCUSSION 



It is interesting to compare results found here for hybrid baryons to the predictions of 
the bag model [yj. Out of all of the flux-tube model states listed under Hf and if-j 4 in 
Table ITVl only the V 2 | + , and V 2 | + states have the same flavor, quark-spin S, total angular 
momentum and parity as the low-lying hybrid baryons in the bag model. However, restrict- 
ing to experimentally measurable quantum numbers (flavor, total angular momentum, and 
parity), only one light-hybrid is different between the flux-tube and bag model predictions. 
This J p = 5/2 + state is flavor A in the flux-tube model and flavor N in the bag model, but 
is amongst the higher-lying states in both models f|. 
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Looking at quantum numbers alone, the Roper resonance could be regarded as a hybrid 
baryon candidate in both models. However, our mass estimates do not support this identifi- 
cation, in contrast to some bag model [f| and QCD sum rule estimates ^|. Both our model 
and the bag model have seven low-lying hybrid baryons. 

One of the disadvantages of hybrid baryons relative to hybrid mesons, multi-quark states 
and glueballs, is that all possible baryon quantum numbers can be attained by conventional 
baryons, so there are no "exotic" quantum numbers. However, our low-lying hybrid 
baryons are "non-relativistic quark model exotic", since these states have quark-label ex- 
change antisymmetric flux wavef unctions, and so totally antisymmetric space x spin x flavor 
wavefunctions, and states of this nature cannot be constructed in the non-relativistic quark 
model. 

The phenomenology of hybrid baryons has been reviewed recently 
be discussed here. 
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331 ] and so will not 



VIII. CONCLUSIONS 

Significant progress has been made towards building a realistic flux-tube model of (hybrid) 
baryons. The full multi-bead Hamiltonian is constructed and it is demonstrated that the 
junction decouples from the beads on the triads to a high degree of accuracy (with the exact 
and decoupled lowest frequencies differing by at most 2%), so that the lowest-lying hybrid- 
baryon excitations in the flux-tube model can be associated with the motion of the junction. 
This simplifies the description of the long distance properties of a low-lying (hybrid) baryon 
to be that of three quarks and a junction, with the junction connected to each of the quarks 
via a linear potential. The parameter dependence of the conventional and hybrid-baryon 
potential is predicted, and can be used as an input in various theoretical approaches. 

The quantum numbers of the lightest hybrid baryons in the flux-tube model are given 
in Table IIV1 and in the presence of the expected spin-spin interactions between the quarks, 
the lightest hybrid baryons are four nucleons with J p = l/2 + , and 3/2 + and with a mass of 
1865± 100 MeV. 
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APPENDIX A: ACTION OF THE PLAQUETTE OPERATOR ON THE JUNC- 
TION 

It is now shown that the junction can move and leave the string in the ground state of its 
Y-shaped configuration in first order HLGT perturbation theory, i.e. with the application 
of a single plaquette operator. 

In Fig.Elthe plaquette operator which operates on the junction has the effect of producing 
a link which is indicated with two arrows on it. The color structure of this link is 3<8>3 = 3©6. 
A link with color 6 is not allowed by conservation of color, because its neighboring link is 
color 3 flowing in the opposite direction, or color 3. However, it would appear that a link 
with color 3 is allowed. The following shows that this is the case. 

An enlargement of the junction region of Fig. El is provided in Fig. ^2 The creation of 
a link from spatial positions A to B with color projections a and b respectively, is denoted 
by Uj^b- Denoting by capital letters the spatial positions, and by lowercase letters the color 
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projections, the mathematical expression corresponding to the graph in Fig. fTTl is 

/] ^ 7 AD L T DG L 1 BG Up 9 £ ^ CF 6 9192S3 ^FE ^ED ^ DG ^GF J (Al) 
dfgig293f'ed'g' 

where the first term is the initial baryon state and the second term is the plaquette operator. 
Now contract the two DG links, and the two FG links, using Appendix B of ref. |34j ] 

U DG U DG = 2 e dd>d>i£ gi gig"U G p + 

d"g" 

T (8dd"8d>d"' + 8dd'"8d'd") (figig"dg'g'" + 8 gig" 1 $g'g") U D q d ^ 99 ^ (A2) 

d"d"'g"g'" 

v'j&v& = ! %«w + 2 £ u °£ v '- < A3 > 

/'V 

where U 6 and £7 8 refer to sextet and octet links respectively, and A a are the usual Gell-Mann 
SU(3) matrices. It follows that Eq. IA1I equals 

-- ^2 U ad u bg U gd U^ f F U F e E UE d D e d d'd", (A4) 

g"fedd'd" 

where the 6 term does not contribute as mentioned above, and the 8 term is not displayed 
as it leads to a more complicated topology. Eq. IA4l is a mathematical expression which gives 
a graph similar to Fig. ^2 with the junction moved from spatial position G to position D, 
and the links in their usual color triplet state. Hence the application of a single plaquette 
operator can move the junction and leave the string in its ground state, as promised. 

APPENDIX B: FLUX-TUBE HAMILTONIAN 

The flux-tube Hamiltonian is first derived for the 6j = 120° case, using the coordinate 
system defined in Eq. |21 The goal is to describe the position of each quark and bead and 
the junction with respect to an alternative coordinate system with origin at the center of 
mass R cm of the entire system. 

In this coordinate system the junction is at rj = r — R cm , where r = (x, y, z). The quark 
positions are 

r t = h(-el,ei,0)-R cm , (Bl) 

where the e l are defined in Eq. Bead n on the z-th triad is at position 

i _ , Nj + l-n ■ n ■ i Nj + l-n { n ■ . n . 

r„-l ^ + 1 e y li+ Ni + 1 x+l; n e x , ^ + 1 ej t + ^ + lV +$ n e v , ^ + ^zj ±t cm . 

(B2) 
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Here ^ is the displacement of bead n along the direction e l which is perpendicular to the 
i-th. triad in the QQQ plane (see Fig. 0), and z l n is its displacement perpendicular to the 
QQQ plane. Note the bead displacements are transverse to the triads on which they lie, and 
are measured from the line connecting quark i to the (in general displaced) junction. Bead 
n = 1 is placed next to the i-th quark, and bead iV, lies next to the junction. 
Requiring that the center of mass is at the origin gives the constraint 

3 Ni 3 

= mjrj + m b J2J2 r n + J2 Mir ^ ( B3 ) 

i=l n=l i=l 

which can be solved for R cm . The Hamiltonian is 

1 1 3 Ni 1 3 3 Ni 

^flux = ^mjPj + -m b + 2 E M ^ + E E ^ ~ <\> ^ 

~ i=l n=l i=l i=l n=0 

where the first three terms form the kinetic energy. The last term is the linear potential 
energy between neighboring beads, where t 1 n+1 is defined as the position of the junction rj 
and Iq is defined as the quark position rj. 

The Hamiltonian is now simplified using the redefined adiabatic approximation, where the 
distances between the quarks remain fixed, concisely stated as U = 0. The small-oscillations 
approximation is used to Taylor expand the potential energy in Eq. IB4I to yield a function 
quadratic in the junction and bead position coordinates. 

The displacements of the beads within and out of the QQQ plane, given by £^ and z l n , is 
expanded in terms of the amplitudes of the m-th normal mode of the beads, q^ m and q]_ m , 
respectively, by 

Ni Ni 

h sr^ i 1717171 i sr^ % 7111171 i r \ 

C = 2^H sm ]vTl Zn = z^ q±mSm W+l' ( ^ 

m=l m=l 

The potential energy now becomes diagonal in these normal mode amplitudes q l ^ m , q]_ m of 
the beads. Carrying out the necessary algebra yields Eq. El As an aside, it is perfectly 
meaningful in Eq. El to put Ni = N 2 = N 3 = 0, the case with no beads other than the 
junction. 

For the 6j > 120° case the definition of the coordinates in Eqs. IB1HB2I remains the same, 
except that the e l are given by Eq. and there are no beads on the third triad. Carrying 
out the necessary algebra yields the Hamiltonian discussed in section III Bl 
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APPENDIX C: QUARK LABEL EXCHANGE SYMMETRY FOR THE 9j = 120° 
CASE IN THE SMALL-OSCILLATIONS APPROXIMATION 



It is now shown, using the explicit definition of r) z} that both quark-label exchange 
symmetric and antisymmetric realizations of rj z are possible. First consider the possibility 
a z = 1 in the definition Eq. "213 of fj z , so that r) z points along p x A. By inserting the 
expressions for p and A in Eq. "2H]into p x A, it can explicitly be verified that p x A is totally 
antisymmetric under exchange symmetry. Hence, f\ z with a z = 1 is totally antisymmetric 
under quark label exchange symmetry. 

Another possibility is to choose fj z along the vector 

Z p x A p x A, (CI) 

where Z is a space-fixed unit vector (not determined by the quark positions). This is 
equivalent to choosing 

= 

\Z-px A| 

in Eq. [213 This choice of a z obviously yields a totally symmetric f] z . 

Next consider the quark-label exchange symmetry of i) and i) in the 9,j = 120° case in 
the small-oscillations approximation, as defined in Eq. "213 Consider the vectors fj' ± , defined 
by f) ± = <J±r)' ± . Applying P^ exchanges Zj <-> lj in Eq. [2H1 as well as the labels i and j in 
x and y in Eq. |2j Under Pi 2 , it is easy to see that fj'± — > — fj' ± . Both Pi 3 and P23 can be 
shown to lead to i)'± — > ±77^., where the sign is dependent on the relative sizes of Zi, Z2, and 
Z3. For example, under P23 one can show that this sign is given by the sign of the expression 

v^iTMaT ±(r + f-f)> ( C3 ) 

'2 '3 f l 

where s(Z l5 Z 2 , Z 3 ) is defined in Eq. 12*71 

The fact that the fj' ± transform under label exchange into themselves, up to a sign, follows 
from the observation that the ray in which the fj'± lie is the physical line of oscillation of the 
junction in these vibrational modes. Since label exchange does not change the physics, the 
oscillation should still be in the same ray after label exchange, as found. Since only the ray 
in which r) ± lies is physical, the possibility cannot be excluded that the fi' ± are multiplied by 
a sign a± when a standard choice of eigenvectors is constructed, as in Eq. "213 It is possible to 



40 



show that a consistent set of sign conventions a± can be chosen such that the fj ± are either 
totally symmetric or antisymmetric under label exchange. Neither choice can be excluded. 
It is also possible to show that the signs cr_, a + and a z are invariant under parity. 
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